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ABSTRACT 


Results are presented for the fully-coupled, radiating flow 
field around the entry probes from a detailed calculation which Includes 
equilibrium chemistry, nongray radiative transport, ablation product 
Injection, and a laminar or turbulent bouadary layer. These results 
show that the radiative flux toward the body Is attenuated In the 
boundary layer at downstream regions of the body as well as at the 
stagnation point and that, even when radiation absorption by ablation 

! 

I 

I products Is accounted for, the radiative heating rates along the down- 

I stream regions of the body can, imder certain conditions, exceed 

i 

the stagnation-point values. It is also shown that, for Venusian 
entry, the spectral distribution of radiative flux and the magnitude 

( 

I of radiation absorption by ablation products depend strongly on 

! entry velocity, and that the state of the boundary layer (l.e., laminar 

1 

I or turbulent) can significantly Influence the amount of ablation pro- 

duct absorption or emission that occurs In various spectral regions. 

! 

i Results from a series of radiative, Invlscid calculations .'re also 

presented. 
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INTRODUCTION 


The proposed multi-probe mission to Venus, Pioneer Venus, is the 
first for which radiative heating is a dominant feature of the aero- 
thermodynamic environment and for which a detailed analysis of radiative 
heating is needed to evaluate heatshleld requirements. The Advanced 
Entry Analysis Branch at the Langley Research Center has a continuing 
program for the development of analytical techniques for calculation 
of the radiating flow field about planetary entry bodies. A recently 
developed technique has been used to obtain several benchmark solutions, 
which are presented herein ^ for conditions near peak heating for the 
Pioneer Venus probes. 

This technique provides for the calculation of the fully-coupled, 
viscous, radiating flow field about an ablating nxisymmetric body. 
Detailed treatments of nongray radiative transport; equilibrium 
chemistry; and mass, momentum and energy transport are Included. The 
flow may be either laminar or turbulent and the ablation rate is 
determined as part of the solution through a steady-state ablation 
calculation. 

Results from a series of Invlscld calculations are also presented. 
These invlscld solutions were obtained in order to determine the 
radiative heating history during entry of the probes, the effect of 
the composition of the Venusian atmosphere on radiative heating, and 
to explore the effect of certain changes in probe geometry on radiative 


heating 
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freest ream enthalpy 
ablation rate of the heat shield 
freestream pressure 
convective heating rate 
radiative heating rate 

momenti::in thickness Reynolds number for boundary layer 

nose radius of body 

base radius of body 

distance along body surface 

trajectory time 

temperature 

entry velocity 

freestream velocity 

weight of entry probe 

ballistic coefficient 

entry angle 

cone angle of body 

freestream density 

aerodynamic shear 


Subscripts : 


e 

w 

o 


refers to the boundary-layer edge 
refers to the body wall 
refers to the stagnation point 
refers to radiation frequency 
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METHOD OF SOLUTION 

A complete description of the method for the fully coupled solu- 
tion of the radiative flow field around an ablating entry body Is 
given in references 1 and 2. Briefly, the method of solution Is to 
separate the flow flald Into an outer layer where the Invlscld flow 
equations are applicable and an Inner viscous layer, either laminar 
or turbulent, where the boundary layer equations are applicable. The 
divergence of the radiative flux Is included In the energy equation 
for the solution of each layer and the boundary layer is coupled to 
the invlscld shock layer by the radiative transport through both 
layers and by the boundary layer displacement thickness. In the 
coupling of the two layers, the Invlscld flow field Is displaced from 
the wall by the boundary layer displacement t.ilckness and the boundary 
layer profiles are used out to the point where the boundary layer 
edge properties and their derivatives equal the corresponding Invlscld 
layer values. The radiative transport for the boundary layer solu- 
tion Is calculated from the matched profiles through the entire 
layer. H'.e treatment of radiation Includes molecular band, continuum, 
and line transitions with a detailed frequency dependence of the 
absorption coefficient. The ablation rate Is computed as part of 
the solution under the assumption of steady-state ablation. 

Several solution options are available. One Is the fully-coupled 
solution described previously. A second Is an Invlscld full flow- 
field solution; and the third an In/lscld stagnation-region solution. 
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RESULTS AND DISCUSSIONS 

Fully-coupled solutions at typical conditions were 
obtained for the entry probes of Pioneer Venus. The entry 
trajectories used In the analysis are pre.>ented in Figure 1 and were 
calculated from nominal entry parameters and body shapes provided by 
the Pioneer Venus Study Team of the Ames Research Center. The atmos- 
pheric model for the Venusian atmosphere Is that presented In refer- 
ence 3 with the gas composition 97 percent carbou dioxide and 3 
percent nitrogen by volume. 

The radiative heating rates at the stagnation point of the 
bodies were calculated along the trajectories to determine the loca- 
tion of peak radiative heating. The "stagnation-point" version of 
the radiating Invlscid flow program was used for these calculations. 

The results are presented In Figure 2 for the large and small probes. 
Even though the nose radius Is smaller, the radiative heating to the 
small probe Is greater due to Its deeper penetration into the denser 
regions of the atmosphere at high velocities. Also, two peaks occur 
In the radiative heating to the small probe. At the higher velocities, 
the thermodynamic state of the gas behind the shock is such that line 
radiation Is a large contributor of the total radiative heating; 
however, the contribution of line radiation decreases as the entry 
probe decelerates and molecular band radiation becomes Important. 

The large probe decelerates at a higher altitude and the radiative 
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^.eating to the probe is almost constant during the time of maximum 
radiative heating. The possibility of a double peak in the time 
history of radiative heating to a probe will depend on the entry 
trajectory. 

Tne flight conditions represented by circular symbols in Figure 
2 were chosen as being representative of peak heating. It is at 
these conditLons for peak radiative heating that fully coupled solu- 
tions with ablation product effects were obtained for the flow around 
the large and small probes. Hereafter, these are referred to as large 
and small probe solutions. The ablation of the heat shield for the 
fully coupled solution is represented by steady-state ablation of a 
phenolic carbon material with elemental mass fractions of carbon, 

0.851; oxygen, 0.110; hydrogen, 0.035; and nitrogen, 0.004. The steady- 
state ablation rates are obtained as part of the fully coupled solu- 
tion and are not prespecified. 

The results from the fully coupled solutions for the large and 
small probes are presented in Figures 3 and 4, respectively. The 
distribution along the body of radiative heating, convective heating, 
ablation rate, aerodynamic shear, and momentum thickness Re 3 rnolds 
number are presented. The spectral distribution of radiative neat 
flux due to continuum processes is also presented for two body loca- 
tions. At these flight conditions, the line radiation is small com- 
pared to the continuum. 


4 
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For the large probe, the ablation rate decreases in the nose 
region and reaches a nearly constant value along the afterbody 
(Fig. 3(d)). As shown in Figure 3(c), the convective heating rate 
decreases along the afterbody but the radiative heating rate Increases 
to values which exceed the stagnation-point value. The radiative 
heating to the wall is greater than the convective heating along the 
entire body. The total heating rate, radiative plus convective, along 
the afterbody is 30 percent to 40 percent less than the stagnation- 
pol \t value. Thus, the decrease in convective heating rate along 
the afterbody counterbalances the effect of the Increase in radiative 
heating rate. The wall temperature resulting from the steady-state 
ablation is 3600° K, within 2 percent, along the entire body. 

For the small probe, there is a decrease in the ablation rate 
along the afterbody (Fig. 4(d)) and both the convective and radiative 
heating rates decrease along the afterbody (Fig. 4(c)). The convective 
heating rate is greater than the radiative heating rate along the 
entire body. This result is opposite to the result obtained for the 
large probe. The total heating rate along the afterbody of the small 
probe is 40 percent to 60 percent less than the stagnation point 
value. The wall temperature resulting from the steady-state ablation 
is 3700° K, within 2 percent along the body. 

The radiative heating rates from Invlscid, radiating solutions 
are coiq>ared with the results for the fully coupled solutions with 


steady^state ablation in Figures 3(b) and 4(b) for the large and 
small probes, respectively. The radiative heating rates from the fully 
coupled solutions are less than the values from the Invlscld, radiating 
solutions by 9 percent and 17 percent, respectively, for the large 
and small probes at the stagnation point of the bodies. These per- 
cent reductions are nearly constant along the bodies. The reduction 
in radiative heating is due to absorption of radiation within the 
boundary layer. 

Spectral distributions of radiat've heating toward the body due 
to continuum processes are presented in Figures 3(g) and 4(g), respec- 
tively, for the large and small probes. (Note, molecular band radia- 
tion is Included with the ' ontlnuum radiation and both are referred 
to as continuum processes) . The spectral distributions are presented 
for the stagnation point and an afterbody location. As shown by the 
results, the absorption of radiation within the boundary layer occurs 
in the spectral range of 5 to 10 eV and is due to self-absorption of 
the fourth positive band system of carbon monoxide, C0(4+). At the 
stagnation point of the small probe there is a slight Increase in 
radiative heating within the boundary layer at 2.5 eV that is due to the 
Swan band system of diatomic carbon, €2(8). A similar effect was 
noted at the other body locations for both the large and small probes, 
but the effect was too small to be presented in the figures. The 
radiative heating due to line processes was less than 10 percent of 
the toted radiative heating at these entry conditions and spectral 
distributions are not presented. 
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The fully coupled solutions assumed laminar flow In the boundary 
layer. Shown In Figures 3(f) and 4(f) are the distributions of 
momentum thickness Reynolds number around the bodies. If the criti- 
cal value for transition to turbulent flow Is about 200, then the bound- 
ary layer would be turbulent along the afterbody for both the large 
and small probes. This could have an appreciable effect on the dis- 
tribution of convective heating, ablation rate, and aerodynamic shear. 

A fully coupled solution for the large probe with a turbulent 
boundary layer was obtained and the results are compared with the 
laminar case In Figure 5. A critical value of 200 for the momentum 
thickness Reynolds number was used for Instantaneous transition to 
turbulent flew. As expected, the convective heating rates, ablation 
rates, and aerodynamic shears Increased with transition to turbulent 
flow. The distribution of radiative heating along the body was 
essentially the same for both laminar and turbulent flow; thus, the 
turbulent boundary layer did not appear to further attenuate the 
radiative heating even though the boundary- layer thickness was as much 
as 80 percent greater than for lam-<nar flow. This phenomenon can be 
explained with the spectral distributions of ...idlatlve flux at the 
wall presented In Figure 5(d) for the turbulent and laminar solutions. 
For the turbulent boundary layer, there Is an Increase In the self- 
absorptlon of the CO (4+) band system In the spectral range of 5 to 8 
eV but this efte.'t Is offset by an Increased emission of the 
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band system at 2.5 eV and the red band system of cyanogen, CN(R), 
at C.6 to 2 eV. The spectral distrl^ution at the boundary-layer 
edge is essentially the same for both the turbulent and laminar boundary 
layers and was presented in Figure 3(g). While the radiative heating 
rates along the afterbody of the large probe are the same for a tur- 
bulent or a laminar boundary layer, the wall spectral distribution 
of the radiation is different for the two cases. It should not be 
inferred that the radiative heating to the body will be the same for 
a laminar or turbulent boundary layer for other entry conditions or 
other body shapes. 

A fully coupled solution was obtained for freestrea iltlons 

representative of a higher entry velocity than that of tL nominal 
trajectories, and the results are presented in Figure 6. The freestream 
conditions and body shape are given in Figure 6(a). The body shape 
is the same as the large probe except the nose radius is slightly 
smaller. Hereafter, this case is referred to as the high velocity 
entry. 

As shown in Figure 6(c), the radiative heating rate is greater 
than the convective heating rate along the entire body. Unlike the 
results previously presented for the large p''*obe, there is only a 
slight increase in the radiative heating along the afterbody. The 
convective heating rate at the stagnation point is slightly less 
than the values just downstream of the stagnation point due 
to an Increased effect of blockage of convective heating by the mass 
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Injection of the ablation products. As shown in Figure 9(d), the 
ablation rate decreases alon^^ the body In the nose region. The 
ablation rate and total heating rate are nearly constant along the 
afterbody. The wall temperature resulting from the steady-scate 
ablation Is 3750° K, within 1 percent, along the body. 

There Is significant absorption of radiation within the boundary 
layer for the high velocity entry as shown b comparison of results 
between the fully coupled solution and an Invlscld, radiating solution 
in Figure 6 (b). The reduction in radiative heating rate is 25 per- 
cent at the stagnation point and 25 percent to 30 percent along the 
afterbody. Speccial distributions of radiative heating toward the 
body due to continuum processes are presente.d In Figure 6 (g) for the 
stagnation point and an afterbody location. For the continuum pro- 
cesses, most of the absorption wlLhln the boundary layer occurs In 
the spectral region of 5 to 10 eV and is due mainly to self-absorp- 
tion of the CO (4+) band. There Is an Increase In radiative flux al: 

2.5 eV and Is due to the € 2 ( 6 ) band At the stagnation point ':;here 
is absorption at photon energies greater than 10 eV and this absorption 
Is attributed to atomic species within the cooler regions of the 
boundary layer. There Is also absorption at the photoionization edge 
of carbon at 8.5 eV. 

The radiative flux due to line processes Is 50 percent of the 
total radiative flux at the nose region of the body and decreases to 
only 20 percent along the afterbody. Spectral distributions of tki 
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radiative heating due to line processes for two locations at the 
nose region are presented In Figure 6(h). Ab.-orptlon of line radiation 
within the boundary layer occurs only for the line groups in the ultra- 
violet at 7.1, 8. A, and 9. A eV. The line radiation in the visible 
and infrared regions of the spectrum, less than 2 eV, is not attenuated 
within the boundary layer. 

Information on the effects of different heatshleld materials on 
the attenuation of radiative heating within the boundary layer vmld 
be useful in the final design of the entry probes. The only heat- 
shieM material considered in the present study was phenolic carbon; 
thus, tlie effect of the heatshleld material on the results cannot be 
determined directly. However, the results at the stagnation point 
for the high velocity case can je conq>ared with prior results pre- 
sented in reference 4 for a phenolic nylon material (elemental mass 
fractions of carbon, 0.61; oxygen, 0.20; hydrogen, 0.10; and nitrogen, 
0,09), The freestream conditions and body nose radius were the same 
as the present high velocity solution except the freestream gas com- 
position was 90 percent CO 2 , 10 percent N 2 for the solution of ref- 
erence 4. A tabulation of the results is presented below. 


Material 

Present 

Phenolic Carbon 

Reference 
Phenolic Ny 

q^, MW/m^ 

26.7 

20.8 

q^, MH/m^ 

36.3 

30.7 

q^, MW/m^ 

3.2 

2.9 

®a» 

0.98 

0.8A 

T , °K 
w* 

3770.0 

3AA0.0 
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The results from the two solutions are in general agreement with 

the ablation rate, convective heating rate, and radiative heating 

rate being slightly lower for the phenolic nylon material. The 

R R 

attenuation of radiative heating within the boundary layer 
2 

was 10 MW/m for both materials. A strong conclusion cannot be 
made about the effect of the heatshleld material due to slight 
differences In computational procedure and freestream gas composition 
for the two solutions. The phenolic nylon material appears to be 
slightly better from the standpoint of lower heating rates and abla- 
tion rate at this entry condition. However, the phenolic nylon may 
be more susceptsble to char failure. 

The exact composition of the Venusian atmosphere is unknown at 
present. The most recent data indicates the composition to be pre- 
dominately CO 2 (up to 100 percent) with the other principal constituent 
being N 2 , The effect of gas composition for CO 2 -N 2 mixtures on the invis- 
cid, radiative heating to the stagnation point of a body is presented in 
figure 7 for two entry conditions. For the higher velocity entry, the radia- 
tive heating rate varies by only 25 percent 0 . er the range of CO 2 content with 
the higher values occurring for the greater percentage of CO 2 . For the 
lower velocity entry, the radiative heating rate can vary by a factor 
of 6 over the range of CO 2 content. This Is due to the domination of 
radiation by the molecular band transitions of C0(4+), CH(V), and 
CN(R) at the lower velocity. For the compositions of present Interest 
for Venusian entry (lees than 10 percent N 2 ) , the effect of composition 
on the stagnation-point, radiative heating rate Is negligible at the 
hlg<. r ' doclty: but, a 22 percent difference occurs with composition 
at the lower velocity. 
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During the evolution of the Pioneer Venus project, two differ- 
ent sets of entry probes have been studied. The first probe set, 
referred to hereafter as Type I probes, consisted of a large probe 
which had a 60° cone angle and three small probes which had 45° cone 
angles. The second set of probes, referred to hereafter as Type II 
probes^ consxsted of a large probe and three small probes, all of 
which had 45° cone angles. In addition to the differences In cone 
angle, there were also differences In nose radii and ballistic coef- 
ficients between the two probe sets. The characteristics of the Type 
I and Type II probes and their respective entry trajectories are 
presented In figure 8. The only small probe considered for each 
probe set Is the one entering at the steepest f iry angle. 

The time histories cf cold-wall, convective heating rate and 
liivlscld, radiative heating for the stagnation point of the probes 
are presented In Figure 9. The cold-wall, convective heating rates 
were calculated by 

q° ^ - 1.88x10“^ P ^ V ^ 

^w,o n '^o> a> 

and this equation Is based on the method presented In reference 5 for 
the heat transfer coefficient. Both the radiative and convective 
heating rates for the large probe are greater for Type II than Type I 
and are due to the deeper penetration into the denser refions of 
the atmosphere at high velocities. For the small probe, the main 
difference Is a decrease in convective heating for Type II and Is 
due to creased nose radius. 
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The Invlscld, radiative heating distributions along the body 
at times near peak heating are presented In Figure 10 for the two 
probe sets. A substantial change in the distribution for the 

large probe resulted from a decrease In conical angle from 60 to 45 
degrees. The radiative heating rate increases along the conical 
region for the 60 ° body but decreases along the conical region for 
the 45° body. While the stagnation point value for the large probe 
Is 60 percent greater for Type II than Type I, the values along the 
flank region (where the heatshleld area is largest) decreased by 
50 percent. 
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CONCLUSIONS 

A conputational method developed for the solution of the 
radiating flow field around planetary entry bodies has been 
used to obtain results for the entry probes of the proposed Pioneer 
Venus mission. 

The method is a detailed calculation for a fully-coupled radiating 
flow field with ablation product Injection into a laminar or turbulent 
boundary layer. The results from the present study for the entry 
probes of Plcneer Venus show that the radiative flux toward the body 
is attenuated In the boundary layer at downstream regions of the body 
as well as at the stagnation point, and that, even when radiation 
absorption by ablation products is accounted for, the radiative heating 
rates along the downstream regions of the body can, under certain 
conditions, exceed the stagnation-point values. It is also shown 
that, for Venusian entry, the spectral distribution vf radiative flux 
and the magnitude of radiation absorption by ablation products depend 
stroi^ly on entry velocity, and that the state of the boundary 
layer (i.e., laminar or turbulent) can significantly Influence the 
amoint of ablation product abeorption or emission that occurs in 
various spectral regions. 

The attenuation of radiation by the boundary layer will reduce 
the radiative heating to the body by up to 20 percent for the 
entry conditions and body shapes which are presently considered 
as nominal for Pioneer Venus. The reduction Increases to 30 percent 


at higher velocity entries 
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Figure 4.- Fully coupled^ radiating flow solution with 
steady-state ablation uf carbon-phenolic 
heat shield for aaall probe. 
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Figure .5.- Comparison between a turbulent and laminar 
boundary layer for fully coupled solution 
for large probe. 
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Figure 7.- Stagnation-point » radiative heating rates In 
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